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Thermal  radiation  in  the  air  and  fuel  channels  has  the  potential  to  dramatically  influence  the  overall 
operating  conditions  and  performance  of  solid  oxide  fuel  cells.  A  3D  comprehensive  model  is  developed, 
with  emphasis  on  quantifying  the  radiative  heat  transfer  process  and  its  effects.  The  radiosity  method 
is  used  for  the  thermal  radiation  in  the  air  and  fuel  channels.  The  thermal  radiation  heat  transfer  is 
coupled  to  the  overall  energy  conservation  equation.  Commercially  available  COMSOL  CFD  software  is 
used  as  a  platform  for  the  global  thermal-fluid  modeling  of  the  SOFC.  The  effects  of  the  operating  voltage, 
emissivity,  ambient  temperature  and  flow  arrangement  (Co-  and  Counter-flow)  on  the  performance  of 
SOFC  are  investigated.  The  predicted  results  reveal  that  the  radiative  heat  transfer  should  be  considered 
in  SOFC  modeling  simulation,  and  the  effects  of  the  thermal  radiation  on  the  performance  of  SOFC  under 
a  different  flow  arrangement  is  sometimes  quite  significant. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  are  promising  energy  sources 
that  are  very  attractive  due  to  their  high  energetic  efficiency  and 
low  emission  of  pollutants  from  the  exhaust  gases  into  the  envi¬ 
ronment.  However,  the  SOFCs  have  their  drawbacks.  There  are 
thermal  expansion  mismatches  among  multifunctional  materials 
due  to  the  high  operating  temperature  (up  to  1200K)  [1].  There¬ 
fore  it  is  of  major  importance  to  predict  accurately  the  temperature 
field  within  the  cells  in  order  to  detect  the  hot  spots  and  the 
maximum  thermal  gradients  [2].  Due  to  the  high  operating  temper¬ 
ature,  radiative  heat  transfer  must  be  given  special  consideration 
in  any  modeling  and  design  efforts,  including  stack  thermal  man¬ 
agement  and  materials  development.  In  SOFCs,  the  radiative  heat 
transfer  usually  consists  of  three  parts:  the  surface-to-surface  heat 
exchange  in  anode  and  cathode  channels;  the  thermal  radiation 
of  participating  media  in  the  reactants  and  products;  the  thermal 
radiation  of  translucent  electrolyte  and  electrode  layers.  During  the 
last  decades,  many  researchers  have  paid  attention  to  investiga¬ 
tions  of  the  radiative  heat  transfer  coupled  with  electro-chemical 
and  transport  processes  in  SOFCs  [3-15]. 

In  references  [3-5],  the  authors  regarded  the  radiative  heat 
transfer  to  have  a  great  effect  on  the  heat  transfer  rate  in  SOFCs. 
Damm  and  Fedorov  [6,7]  investigated  the  effects  of  radiative  heat 
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transfer  in  the  electrode  and  the  electrolyte  layers.  Their  results 
revealed  that  the  electrodes  can  be  regarded  as  optical  transpar¬ 
ent  material,  and  the  radiative  heat  transfer  within  the  electrodes 
had  a  negligible  effect  on  the  average  cell  operating  temperature, 
voltage,  or  temperature  gradients.  However,  the  effects  of  radia¬ 
tion  heat  transfer  within  the  electrolyte  depend  on  the  thickness  of 
electrolyte  layer,  i.e.,  the  thicker  the  electrolyte  layer,  the  greater 
the  impact  of  radiative  heat  transfer.  The  radiative  heat  transfer 
had  little  influence  when  the  thickness  of  the  electrolyte  was  less 
than  1 5  pum.  The  same  conclusion  was  drawn  by  Daun  et  al.  [8].  The 
radiative  heat  transfer  with  the  participating  media  in  SOFCs  was 
investigated  in  detail  by  Vandersteen  and  Pharoah  [9].  Because  of 
the  low  emissivity  of  the  gas  and  the  small  channel  sizes,  the  gas 
was  optically  thin,  and  the  effect  of  participating  media  on  thermal 
radiation  was  minimal  in  the  planar  geometry,  but  it  was  likely  to 
be  significant  in  the  tubular  geometry.  Sanchez  et  al.  [  1 0]  also  inves¬ 
tigated  the  radiative  heat  transfer  with  participating  media  in  the 
anode  and  they  found  it  had  little  influence  on  the  performance  of 
SOFCs  under  the  ordinary  operation  conditions.  They  also  indicated 
that  the  surface-to-surface  thermal  radiation  in  the  flow  channels 
had  great  effects  on  the  SOFCs.  Yakabe  et  al.  [1 1  ]  numerically  inves¬ 
tigated  the  effects  of  surface-to-surface  thermal  radiation  on  the 
performance  of  planar  SOFCs,  and  revealed  that  the  temperature 
distribution  in  the  cell  became  flat  as  the  radiative  heat  exchange 
inside  the  channels  was  considered. 

The  radiative  heat  transfer  has  received  much  more  attention  in 
recent  years  [2,12-14].  Based  on  a  Monte  Carlo  Ray  Tracing  method, 
a  numerical  methodology  aimed  to  predict  the  thermal  radiation  of 
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Nomenclature 

Av  relative  surface  area  per  unit  volume  (m2  m-3 ) 
c  concentration  of  species  (mole  m-3 ) 

cp  specific  heat  (J  kg-1  K-1 ) 

d  pore  diameter  (m) 

effective  diffusion  coefficient  (m2  s) 

F  view  factor 

G  irradiation  ( W  m-2 ) 

h  convective  heat  transfer  coefficient  (W  m-2  K-1 ) 

H  channel  height  or  layer  thickness  (m) 

i  current  density  (A  nrr 2 ) 

J  radiosity  (Wm-2);  kinetically  controlled  local  cur¬ 

rent  density  (A  m-3 ) 

L  cell  length  (m) 

m  molar  flow  rate  (mol  s-1) 

n  number  fraction  of  particles 

p  pressure  (Pa);  probability  of  percolation  in  reaction 

zone 

qr  net  heat  transfer  rate  from  a  surface  ( W  nrr 2 ) 

r  radius  of  electron-conduction  particles  (m) 

T  temperature  (I<) 

u  velocity  vector  (ms-1) 

W  flow  channel  width  (m) 

x  molar  fraction 

Greek 

s  porosity;emissivity 

iq  overpotential  (V) 

0  contact  angle  between  electron-  and  ion-  conduct¬ 

ing  particles  in  reaction  zone  (°) 
k  permeability  (m2) 

A  thermal  conductivity  (W  m-1 1<-1 ) 

a  electrical/ionic  conductivity 

(S  m-1 );  Stefan-Boltzmann  constant 
(5.67  x  10-8  Wm-2  K-4) 
r  tortuosity 

cp  potential  (V) 

0  volume  fraction 

(Di  mass  fraction  of  species  i 

Subscripts 

amb  ambient 

ASL  anode  support  layer 

ARL  anode  reaction  layer 

CRL  cathode  reaction  layer 

E  electrolyte 

CDL  cathode  diffusion  layer 

el  electronic 

io  ionic 

Superscrits 

eff  effective 


Table  1 

Geometry  of  unit  cell. 


Geometry 

Value 

Reference 

Cell  length!  (m) 

2.0  x  10“2 

/ 

Flow  channel  width  WCh  (m) 

1.0  x  10“3 

[17] 

Flow  channel  height  Hc h  (m) 

1.0  x  10“3 

[17] 

Anode  support  layer  thickness  HA SL  (m) 

1.0  xl0“3 

[16,17] 

Anode  reaction  layer  thickness  HArl  (m) 

2.0  x  10“5 

[16,17] 

Electrolyte  thickness  HE  (m) 

8.0  xlO-6 

[16,17] 

Cathode  reaction  layer  thickness  Hc rl  (m) 

2.5  x  10“5 

[16,17] 

Cathode  diffusion  layer  thickness  Hc dl  (m) 

1.3  x  10“5 

[16,17] 

modeling  was  introduced  by  Bao  et  al.  [14].  The  detailed  radiation 
model  based  on  analytical  view  factors  predicted  more  uniform 
distribution  of  the  cell  temperature  and  current  density  in  the 
overall  SOFC  modeling.  Nevertheless,  most  of  the  existing  heat 
transfer  models  for  the  SOFCs  simply  ignored  the  effect  of  thermal 
radiation  [15],  even  though  this  effect  is  important  when  the  cell 
operating  temperature  is  higher  than  800  K  [4].  As  some  authors 
argued,  the  view  factor,  depending  on  the  cell  geometry  and  con¬ 
figuration,  has  a  great  influence  on  the  surface-to-surface  radiation 
exchange.  However,  2D  models,  available  in  most  existing  studies 
[4,8,10,12,15],  could  not  account  for  this  phenomenon. 

By  considering  various  configurations  of  interconnect,  fuel 
channel,  anode  support  layer,  anode  reaction  layer,  electrolyte, 
cathode  reaction  layer,  cathode  diffusion  layer  and  air  channel,  a  3D 
fuel  cell  unit  model  is  developed,  and  implemented  in  the  software 
COMSOL  to  verify  if  radiative  heat  transfer  is  really  an  important 
phenomenon  in  flow  channels  of  the  typical  planar  SOFCs.  Next,  we 
validate  the  model  with  some  experimental  data.  Finally,  with  the 
validated  model,  we  investigate  the  effects  of  emissivity,  ambient 
temperature  and  flow  configuration  on  the  SOFC  performance. 

2.  Physical  model  and  computational  method 

2.1.  Physical  model 

In  general,  one  cell-stack  is  constructed  consisting  of  repeated 
single  unit  cells.  Fig.  1(a)  shows  a  schematic  diagram  of  single 
cell  composed  of  five  pairs  of  parallel  gas  flow  channels.  Since 
much  more  time  is  needed  to  simulate  the  cell  performance  for 
the  whole  stack,  a  simplified  single-unit  cell  with  all  the  compo¬ 
nents  is  selected  in  the  present  work.  Fig.  1(b)  and  (c)  shows  the 
repeated  unit  cell  located  in  the  middle  part  of  the  stack.  For  vali¬ 
dation  of  the  present  model,  the  configuration  and  the  parameters 
are  chosen  based  on  the  experimental  work  by  Jung  et  al.  [16].  The 
detailed  geometry,  operating  parameters  and  properties  are  listed 
in  Tables  1-3,  separately. 

2.2.  Mathematical  formulations  and  computational  method 
2.2.1.  Thermal  radiation  model 

Inclusion  of  radiative  heat  transfer  in  the  analysis  of  the  cell 
performance  includes  a  number  of  challenges,  particularly  as  the 


materials  used  in  the  cell  design  with  a  planar  geometry  was  devel¬ 
oped  in  [2].  Kee  et  al.  [12]  developed  a  relatively  simple  model  to 
rapidly  evaluate  various  configurations  and  operating  conditions 
for  tubular  anode-supported  SOFC  stacks,  and  they  concluded  that 
the  radiative  heat  transfer  is  remarkably  effective  at  removing  the 
heat  from  tube  bundles  of  the  stack.  A  surface-to-surface  radiation 
model  was  employed  in  [13]  to  analyse  the  influence  of  different 
operating  conditions  on  the  temperature  distribution  in  the  anode- 
supported  tubular  cell.  A  mathematical  evaluation  of  view  factors 
for  radiative  heat  exchange  in  longitudinally  distributed  SOFC 


Table  2 

Baseline  conditions. 


Variable 

Value 

Reference 

Operating  temperature  T  (I<) 

1073 

[16,17] 

Operating  pressure  p  (Pa) 

1.013  x  105 

[16,17] 

Anode  gas  inlet  molar  flow  rate  mA\  (mol  s-1 ) 

5.68  x  10“5 

[17] 

Cathode  gas  inlet  molar  flow  rate  mc i  (mol  s-1 ) 

1.48  xlO”4 

[17] 

Inlet  molar  fraction  of  H2  at  anode  *h2i 

0.97 

[16,17] 

Inlet  molar  fraction  of  H20  at  anode  Xh2oi 

0.03 

[16,17] 

Inlet  molar  fraction  of  02  at  cathode  Xo2i 

0.21 

[16,17] 

Inlet  molar  fraction  of  N2  at  cathode  *n2i 

0.79 

[16,17] 
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(  a  )  3D  schematic  diagram  of  cell-stack  (  b  )  Schematic  diagram  of  yz  plane  and  boundary  nodes 


(  c  )  Schematic  diagram  of  xz  plane 


Fig.  1.  Computational  domain. 


thermal  radiation  is  coupled  with  analysis  of  convective- 
conductive  heat  transfer  [7,21  ].  The  first  is  the  inherent  complexity 
of  the  radiation  governing  equations,  which  are  integral- 
differential  equations  and,  in  general,  depend  on  as  many  as 
seven  independent  variables  (time,  three  position  variables,  two 
angular  variables  describing  direction  of  propagation  of  radiation 
rays,  and  the  wavelength).  Furthermore,  the  governing  equations 
are  non-linear,  as  the  emissive  power  features  a  fourth-power 
dependence  on  the  temperature.  Besides  the  difficulty  associated 
with  solving  these  equations,  the  accuracy  of  any  analysis  is  always 
limited  by  the  extent  to  which  radiative  properties  are  known 
[7].  Many  researchers  tried  to  use  different  methods  to  model  the 
radiation  heat  transfer,  such  as  the  Monte  Carlo  method  [2,8,9], 
Schuster-Schwartchild  two-flux  approximation  method  [4,6,8], 
and  Hottel  zone  method  [5].  To  make  a  simple  radiation  analysis 
possible,  all  the  surfaces  in  the  channels  are  assumed  to  be  opaque, 
diffuse  and  gray  in  this  study.  The  radiosity  method  is  introduced  to 
determine  the  net  radiation  heat  transfer  rate  of  the  surfaces.  This 
method  includes  the  concepts  of  radiostiy  and  irradiation,  etc. 

The  radiosity  represents  the  rate  at  which  radiation  energy 
leaves  a  unit  area  of  a  surface  in  all  directions  [21  ].  For  an  opaque 
surface,  the  radiosity  can  be  expressed  as 

J  =  (1  -  e)C  +  ScrT4  (1) 

where  J  is  the  radiosity  (Wm-2);  s  is  the  surface  emissivity;  G 
is  the  irradiation  (Wm-2);  o  is  the  Stefan-Boltzmann  constant 
(5.67  x  10-8  W m-2  I<-4). 

The  radiation  flux  incident  on  a  surface  from  all  directions  is 
called  irradiation  G,  and  is  expressed  as  [21  ] 

G  =  Cm+FambcrT4amb  (2) 


where  Gm  is  the  mutual  irradiation,  coming  from  other  surfaces  in 
the  model  (Wm-2);  Famb  is  the  ambient  view  factor;  Tamb  is  the 
ambient  temperature  (I<). 

For  the  specific  point,  m(x,  y,  z),  in  a  surface,  the  irradiation  of 
the  point  m’  {x\  y' ,  z')  to  the  point  m(x,  y,  z)  can  be  defined  as 


Gm  — 


[  (~n'  •  r)(n  •  i 

Js  W4 


-yds 


(3) 


where  S  is  the  wall  surface  area  (m2);  n'  is  the  unit  normal  factor 
of  m'  (x',  y\  z')  in  solid  surface,  r  is  the  vector  from  the  point  m(x, 
y,  z)  to  m'  (x',  y',  z’);  n  is  the  unit  normal  factor  of  m(x,  y,  z)  in  solid 
surface;  /  is  the  radiosity  of  m'  (x',y',  z'). 

In  Eq.  (2),  the  ambient  view  factor  is  evaluated  by 


Famb  =  l—  /  (~n'  'r)("  r)dS  (4) 

amb  Js,  tt \r\4 

According  to  the  energy  balance,  the  net  flux  of  the  heat  transfer 
from  a  surface  is  denoted  by  qr  and  is  expressed  as 


Qr  =J  -  G 


(5) 


With  Eqs.  (1 ),  (2)  and  (5),  the  net  radiation  heat  transfer  rate  can 
be  acquired. 


2.2.2.  Governing  equations  and  boundary  conditions 

Transport  phenomena  in  SOFCs  contain  momentum,  mass,  elec¬ 
tron,  ion  and  heat  transfer,  which  are  coupled  with  each  other. 
Due  to  the  multiphysics  nature  of  SOFCs,  various  governing  equa¬ 
tions  (except  for  continuity  equation)  and  boundary  conditions  are 
outlined  separately  in  the  following  sections. 
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Table  3 

Properties  of  material. 


Variable  Value  Reference 


Porosity  of  anode  support  layer  £Asl 

0.42 

[16,17] 

Porosity  of  anode  reaction  layer  £ARL 

0.097 

[16,17] 

Porosity  of  cathode  support  layer  £Crl 

0.40 

[17] 

Porosity  of  cathode  diffusion  layer  £Cdl 

0.36 

[16,17] 

Average  pore  diameter  of  anode  support  layer  dASL  (pm) 

1.40 

[16,17] 

Average  pore  diameter  of  anode  reaction  layer  dARL/(pm) 

0.33 

[16,17] 

Average  pore  diameter  of  cathode  reaction  layer  dcRL  (pm) 

2.00 

[16,17] 

Average  pore  diameter  of  cathode  diffusion  layer  dCDL(pm) 

1.40 

[16,17] 

Tortuosity  of  anode  support  layer  rt0rt,ASL 

3 

[17] 

Tortuosity  of  anode  reaction  layer  rt0rt,ARL 

3 

[17] 

Tortuosity  of  cathode  reaction  layer  rt0rt,CRL 

3 

[17] 

Tortuosity  of  cathode  diffusion  layer  rt0rt,cDL 

3 

[17] 

Permeability  of  anode  support  layer  kAs l  (m2) 

3.4  x  10-14 

[16,17] 

Permeability  of  anode  reaction  layer  /cARL  (m2) 

3.4  x  10“14 

[17] 

Permeability  of  cathode  reaction  layer  kqrl  (m2) 

5.4  x  10“14 

[16,17] 

Permeability  of  cathode  diffusion  layer  /ccm/m2 

3.7  x  10-14 

[16,17] 

Volume  fraction  of  electron-conducting  particles  (Ni)  in  anode  support  layer  d>eiASL 

0.4 

[17] 

Volume  fraction  of  electron-conducting  particles  (Ni)  in  anode  reaction  layer  <£>eiARL 

0.4 

[17] 

Volume  fraction  of  electron-conducting  particles  (LSM-YSZ)  in  cathode  reaction  layer  d>eicRL 

0.587 

[17] 

Volume  fraction  of  electron-conducting  particles  (LSM)  in  cathode  diffusion  layer  d>elCDL 

1.0 

[17] 

Electrical  conductivity  of  anode  support  layer  crelASL  (S  irr 1 ) 

1.011  x  105 

[16,17] 

Electrical  conductivity  of  anode  reaction  layer  cr elARL(S  m- 1 ) 

1.011  x  105 

[16,17] 

Ionic  conductivity  of  anode  reaction  layer  <rioARL  (S  m-1 ) 

(1  -  <PARl)  | 

(  1  -  £arl  ^  _ 

CTioEL 

\  Xtort,ARL  J 

/-I  0300  \ 

[18] 

Ionic  conductivity  of  electrolyte  layer  crioE  (S  irr 1 ) 

3.34  x  10  exp  1  — - —  1 

[19] 

Ionic  conductivity  of  cathode  reaction  layer  crioCRL  (S  m-1 ) 

(1  -  @CRl)  | 

(  1  -  £crl  ^  _ 

i  1  ^ioEL 

\  Xtort.CRL  J 

[18] 

Electrical  conductivity  of  cathode  reaction  layer  crelCRL  (S  m_1 ) 

9.3  xlO3 

[16,17] 

Electrical  conductivity  of  cathode  diffusion  layer  <TeicDL(S  m-1 ) 

1.52  x  104 

[16,17] 

Electrical  conductivity  of  interconnect  creiiN  (S  m- 1 ) 

1.13  x  106 

/ 

Contact  angle  between  electron-  and  ion-  conducting  particles  in  reaction  layer  0(°) 

30 

[17] 

Radius  of  electron-conduction  particles  rd  (m) 

i-'- 

o 

X 

o 

[20] 

Radius  of  ion-conduction  particles  rio  (m) 

i'' 

o 

X 

o 

[20] 

Thermal  conductivity  of  anode  kA  (W  irr 1  K-1 ) 

6 

[20] 

Thermal  conductivity  of  electrolyte  kE  (W  m-1  K-1 ) 

2.7 

[17] 

Thermal  conductivity  of  cathode  kc  (W  irr 1  K_1 ) 

11 

[17] 

(1)  Momentum  equations 

In  the  air  and  fuel  channels, 


vector  (m/s),  p  the  pressure  (Pa).  De^  is  an  effective  diffusion 
coefficient,  which  is  calculated  as  follows  [22,23], 


-V  •  [/*(Vu  +  (Vu)7)]  +  p{ U  ■  V)u  +  Vp  =  0  (6) 

In  the  porous  electrodes,  the  Brinkman  equation  is  adopted, 

(^+Sm)u=iv.[At(Vu  +  (Vu)7)]-Vp  (7) 

where  /x  is  the  dynamic  viscosity  of  the  gas  mixture;  /c  the  per¬ 
meability  of  the  porous  media;  £  the  porosity  of  the  porous 
media;  Sm  the  source  term  of  continuity. 

The  inlet  condition  for  air  and  fuel  is  as  follows 


u  =  Uq  •  n 


(8) 


Df 


J_\_1 

^  \  DjK  J 


(11) 


where  Dy  is  the  binary  diffusion  coefficient  and  DiK  denotes 
Knudsen  coefficient  for  the  species  i. 

The  source  term  Sz  denotes  the  production  or  consumption 
of  the  species  i  due  to  electrochemical  reactions,  which  can  be 
expressed  as: 

In  the  cathode  reaction  layer, 


(12) 


In  the  anode  reaction  layer, 


The  outlet  boundary  condition  for  air  and  fuel  is  set  as  a  pres¬ 
sure  outlet. 


Sh2 


JMu2  ~  _  JM h2o 

2 F  ’  H2°  2 F 


(13) 


P  —  Po  (9) 

(2)  Mass  transfer  equation 

The  Maxwell-Stefan  equation  involving  Knudsen  diffusion  is 
used, 

V  (-po,-  J2 Df  (Vx*  +  (*>  -  ®i)-y )  +  wu)  =  s>  (1°) 

where  p  denotes  the  density  (kg/m3),  coi  mass  fraction  of 
species  t,  Xj  the  molar  fraction  of  species  j,  u  the  velocity 


The  inlet  mass  fraction  for  air  and  fuel  is  specified. 

cot  =  co  uo  (14) 

At  the  outlet,  the  convection  is  the  dominating  process  to 
force  the  mass  flow  through  the  channel  exits,  and  the  mass  flux 
due  to  the  mass  diffusion  across  the  outlet  is  zero.  The  following 
condition  is  applied  to  all  species. 

n  ■  (- pcoi  J2  Df  (Vx<  +  (*>  -  )  =  0  (1 5) 

(3)  Charge  transfer  equations 
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Conservation  of  electronic  charge  transfer  reads, 
V-(-oreiV^el)  =  Sel  (16) 

Conservation  of  ionic  charge  transfer  is, 

V-(-aioV^io)  =  Sio  (17) 

In  the  cathode  reaction  layer, 

Sel  =  -Si0=J  (18) 

In  the  anode  reaction  layer, 

Sei  =  -Si0  =  -J  (19) 

In  the  electrolyte  layer, 

Si0  =  0  (20) 

The  boundary  conditions  for  the  charge  transfer  are  outlined 
as  below.  As  shown  in  Fig.  1(b),  the  electric  ground  for  Eq.  (16) 
is  set  at  the  top  wall, 

<pel  =  0  (21) 

For  boundary  30  shown  in  Fig.  1  (b),  the  operating  voltage  Vcen 
is  specified 

(pe\  =  Vceu  (22) 


The  remaining  boundaries  for  electronic  and  ionic  transport 
are  set  as  an  insulation  condition, 


where  0  is  the  contact  angle  between  the  electron-  and  ion¬ 
conducting  (Ni  and  YSZ)  particles  in  the  reaction  zone  layers, 
rel  the  radius  of  the  electron-conducting  (Ni)  particles,  nt  the 
total  number  of  particles  per  unit  volume,  nel  and  nio  are  the 
number  fractions  of  the  electron-  and  ion-conducting  particles 
in  the  reaction  zone  layers,  respectively.  Zel  andZio  are  the  coor¬ 
dination  numbers  of  the  electron-  and  ion-conducting  particles 
in  the  reaction  zone  layers,  respectively,  Z  is  the  total  average 
number  of  contacts  of  each  particle,  pei  andpio  are  the  probabili¬ 
ties  of  the  electron-  and  ion-conducting  particles  in  the  reaction 
zone  layers,  respectively.  The  detailed  discussion  on  the  above 
mentioned  parameters  can  be  found  in,  e.g.,  Refs.  [25,26,18]. 

(5)  Energy  equation 

The  local  thermal  equilibrium  assumption  is  adopted  in  the 
heat  transfer  model. 

V  •  (-A.VT)  =  Sq-  pCpU  •  VT  (29) 

The  equation  is  applied  for  the  whole  domain,  but  the  convec¬ 
tion  term  is  considered  only  for  the  fuel  and  oxidant  channels, 
and  omitted  in  the  porous  electrode  for  the  velocity  in  the  elec¬ 
trodes  is  quite  small.  The  source  termSq  in  Eq.  (29)  is  formulated 
as: 

in  the  interconnect  and  anode/cathode  diffusion  layer, 

Sq  =  <?e\V<Pe\  '  ^<Pe I  (30) 

in  the  electrolyte, 


n  •  (-crV<p)  =  0 


(23)  Sq  =  o-ioV<pio  ■  V<pi0 


(31) 


(4)  Electrochemistry  equations 

The  source  term  for  Eqs.  (16)  and  (17)  can  be  acquired  by 
Butler-Volmer  equations  as 


Jan  —  AJq \ 


ref 


_  \  Yh2 
ch2  \ 

Ch2,  ref  J 


X  {exp  (anFy-an  )  -  exp  (-(1  ~  )  }  (24) 


x  {exp  (ggfltog)  _  exp  (J1  -  “gW)  }  (25) 


where  Av  is  the  actual  reactive  surface  area  per  unit  volume, 
Jq2  and  are  the  reference  exchange  current  densities 
for  H2  oxidation  and  02  reduction  reactions  at  the  reference 
concentrations,  cH2?re/  and  c02,re/,  respectively,  and  a  is  the 
charge  transfer  coefficient  (or  symmetry  factor),  whose  value 
lies  between  zero  and  unity,  n  the  number  of  electrons  partic¬ 
ipating  in  the  electrochemical  reactions  and  r)act  the  electrode 
activation  overpotential  (anode  or  cathode),  defined  as  [24] 

ljact,an  =  Vrev,an  ~  |<Pel  —  ^io  |  (26) 

dact,ca  =  Vrev,ca  ~  |<PeI  —  ^io  |  (27) 

As  stated  in  Ref.  [25],  the  values  of  cH2,re/  and  c02,re/  are 
10.78  mol  m-3  and  2.38  mol  m-3,  and  the  corresponding  val¬ 
ues  of  the  reference  exchange  current  densities  Jj1  re/  anc* ■/? re/’ 
are  1320  (Am-2)  and  400  (Am-2),  respectively.  In  order  to 
enhance  the  predictive  capability  of  the  model,  the  expression 
used  to  model  the  reactive  surface  area  per  unit  volume  is  based 
on  random  packing  of  binary  spherical  particles  developed  by 
Costamagna  et  al.  [26].  This  is  given  as 


Av  =  tt  sin2  0r^ntneinio 


Ze\Zjp 

Z 


PelPio 


(28) 


in  the  anode/cathode  reaction  layer, 

Sq  =  adV<pd  ■  V(pd  +  (Ti0Vipi0  •  V<pi0  -  jpTAS -]r)act  (32) 
For  the  inlet,  the  temperature  of  air  and  fuel  is  given  as 


T  =  T0  (33) 

For  the  outlet,  only  the  convective  heat  tranfer  is  considered 
and  local  one-way  flow  is  assumed. 

q  =  n  (pCpUT)  (34) 

For  the  air  and  fuel  channel  walls,  the  net  radiation  heat  flux 
qx  is  also  included. 


n.(-A.VT)  =  /i(Tw-T/)  +  qr  (35) 

where  the  convective  heat  transfer  coefficient  h  is  deduced  from 
the  Nusselt  number  which  is  recommended  as  3.09  in  Ref.  [27]. 


NllXgas 


(36) 


The  hydraulic  diameter  is  defined  as, 


2WChH(Ih 
WCH  +  HCh 


(37) 


2.2.3.  Model  validation 

Before  proceeding  further,  it  is  necessary  to  ascertain  the  reli¬ 
ability  and  accuracy  of  the  present  numerical  model  and  the  code 
developed.  Grid  independence  tests  are  carried  out  under  the  con¬ 
dition  of  operating  voltage  0.75  V,  emissivity  0.5,  and  the  cell 
length  0.02  m.  A  tetrahedral  mesh  is  adopted.  Fig.  2  shows  the 
variation  of  the  average  outlet  temperature  of  air  with  different 
numbers  of  grid  points.  When  the  grid  mesh  number  changes 
from  38,700  to  103,200,  the  average  outlet  temperature  of  air 
changes  from  11 08.8  K  to  11 10.4 1<  (i.e.,  0.14%).  With  the  increase 
of  the  number  of  points,  the  increase  rate  of  average  outlet  tem¬ 
perature  decreases,  and  the  average  outlet  temperature  varies 
slightly  when  the  number  of  grid  points  reaches  103,200,  which 
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Mesh 

Fig.  2.  Grid  independence  test. 


is  used  for  taking  both  the  accuracy  and  convergence  rate  into 
account. 

In  order  to  validate  the  present  model,  quantitative  compar¬ 
isons  are  also  made  with  the  cases  investigated  experimentally  in 
[16],  where  the  cell  length  is  5  cm.  The  comparison  results  with 
and  without  thermal  radiation  are  shown  in  Fig.  3.  In  the  radiation 
model,  the  emissivity  is  0.5.  When  the  current  density  is  among 
0-3000  Am-2,  the  predictions  by  both  models  have  a  certain  devi¬ 
ation  from  the  experimental  results.  When  the  current  density  is 
greater  than  3000  Am-2,  the  numerical  result  with  the  radiation 
model  agrees  well  with  the  experimental  data.  When  the  current 
density  is  greater  than  15,000  Am-2,  the  simulation  result  with¬ 
out  radiation  model  deviates  from  the  experimental  data  gradually. 
It  is  clear  that  the  current  density  predicted  by  the  model  with¬ 
out  radiation  is  smaller  than  the  one  with  the  radiation  at  the 
same  voltage.  This  means  that  the  model  without  thermal  radiation 
underestimates  the  performance  of  SOFCs.  This  is  so  because  with 
the  increase  of  the  current  density,  the  cell  average  temperature 
and  the  temperature  difference  between  the  inlet  and  the  outlet 
increase  quickly,  and  radiative  heat  transfer  becomes  increasingly 
significant,  which  leads  to  the  difference  in  the  calculation  results. 
The  comparison  result  illuminates  that  present  model  with  thermal 
radiation  is  convincing. 


3.  Results  and  discussion 

3.  I .  Influence  of  emissivity  on  SOFC  performance 

The  emissivity  of  a  surface  represents  the  ratio  of  the  radiation 
emitted  by  the  real  surface  to  that  emitted  by  a  blackbody  at  the 
same  temperature.  Emissivity  is  a  measure  of  how  closely  a  real 
surface  approximates  a  blackbody.  The  emissivity  of  a  real  surface 
is  not  a  constant.  Rather,  it  varies  with  the  temperature  of  the  sur¬ 
face  as  well  as  the  wavelength  and  the  direction  of  the  radiation 
[21  ].  A  brief  survey  of  the  literature  reveals  that  the  thermal  radia¬ 
tive  properties  of  the  materials  used  in  a  SOFC  are  not  easy  to  find 
at  their  operating  temperatures  [2],  and  it  is  difficult  to  describe 
the  emissivity  of  the  channel  and  electrode  surfaces  exactly.  In  the 
present  study,  the  same  emissivity  of  the  channel  and  electrode 
surfaces  is  employed  similar  to  what  was  adopted  in  [4]. 

Figs.  4  and  5  show  the  local  temperature  distribution  for  the 
cross-section  at  the  positions  close  to  the  inlet  and  outlet,  in  which 
the  maximum  and  minimum  temperatures  are  indicated.  Table  4 
lists  the  average  wall  temperature  in  the  air  channel  and  the  max¬ 
imum  and  minimum  temperatures  for  the  cross-section  at  the 
positions  close  to  the  inlet  and  outlet.  In  Fig.  4  and  Table  4,  it  is  found 


Table  4 

Effect  of  emissivity  on  the  average  temperature  in  different  x  plane  (Unit:  K). 


e 

Temperature  at  top  and  bottom  air  channel  walls 

The  maximum  and  minimum  temperature  in  cross  section 

T/vall-23 

Fwall-27 

AT  -  Twaii 23  -  Twau 27 

Tmin 

Tmax 

AT-  Tmax  —  Tm  in 

0.00 

1086.60 

1086.00 

0.60 

1085.36 

1086.71 

1.35 

0.05 

1085.26 

1084.66 

0.60 

1084.11 

1085.38 

1.27 

0.10 

1085.01 

1084.41 

0.60 

1083.88 

1085.12 

1.24 

_  0.20 

1084.75 

1084.16 

0.59 

1083.64 

1084.87 

1.23 

x  =  0.002  m  04Q 

1084.51 

1083.92 

0.59 

1083.42 

1084.63 

1.21 

0.60 

1084.38 

1083.79 

0.59 

1083.30 

1084.49 

1.19 

0.80 

1084.28 

1083.71 

0.57 

1083.21 

1084.40 

1.19 

1.00 

1084.21 

1083.65 

0.56 

1083.15 

1084.32 

1.18 

0.00 

1124.06 

1123.42 

0.64 

1122.95 

1124.19 

1.24 

0.05 

1116.63 

1116.02 

0.61 

1115.57 

1116.76 

1.19 

0.10 

1114.89 

1114.29 

0.60 

1113.84 

1115.01 

1.17 

0.20 

1113.18 

1112.59 

0.59 

1112.14 

1113.30 

1.16 

x  =  0.01 8  m  040 

1111.57 

1110.99 

0.58 

1110.53 

1111.69 

1.16 

0.60 

1110.69 

1110.13 

0.56 

1109.66 

1110.81 

1.15 

0.80 

1110.11 

1109.55 

0.56 

1109.08 

1110.23 

1.15 

1.00 

1109.68 

1109.13 

0.55 

1108.65 

1109.79 

1.14 
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(a)  £=0  ,(b)  £=0.2,(c)  £=0.6,(d)  £=1.0 

Fig.  4.  Local  temperature  distribution  in  the  plane  of  x  =  0.002  m  (close  to  the  intlet,  Co-flow,  unit:  K). 


that  the  temperature  distribution  profiles  are  almost  the  same  for 
the  cases  without  and  with  radiation,  e.g.,  the  maximum  temper¬ 
ature  appears  near  the  electrolyte  and  the  minimum  one  in  the 
middle  region  of  the  air  channel.  With  the  increase  of  emissivity,  the 
maximum  and  minimum  temperatures  decrease  slightly,  as  shown 
in  Fig.  4(b-d),  and  the  temperature  difference  between  the  maxi¬ 
mum  and  minimum  temperature  decreases  as  well  (Table  4).  On 
the  other  hand,  the  temperaure  difference  between  the  top  wall 
and  the  bottom  wall  in  the  cathode  is  almost  the  same.  As  for  the 
postion  close  to  the  outlet  shown  in  Fig.  5  and  Table  4,  the  lowest 
temperature  appears  in  the  anode  side  instead,  and  both  the  max¬ 
imum  and  minimum  temperatures  are  higher.  It  is  also  clear  that 
there  is  a  great  difference  in  temperatures  for  the  cases  without 
(Fig.  5(a))  and  with  radiation  (Figs.  5(b-d)).  Flowever,  the  temper¬ 
ature  difference  between  the  walls  is  roughly  the  same  (0.6  K),  and 
the  variation  of  the  difference  between  the  highest  and  the  lowest 
temperatures  is  kept  almost  similar  for  the  specific  cases  studied. 
Fig.  6  shows  the  local  temperature  distribution  at  the  same  oper¬ 
ating  voltage  (Vcen  =  0.75  V).  The  trends  of  the  local  temperature 
distribution  for  all  the  cases  are  similar,  i.e.,  the  local  temperature 
increases  along  the  main  flow  direction  due  to  the  exothermic  elec¬ 
trochemical  reaction.  Flowever,  with  the  increase  of  emissivity,  the 
local  temperature  decreases,  which  leads  to  a  decrease  in  the  local 
current  density  compared  to  the  case  without  radiation,  as  shown 
in  Fig.  7.  When  the  temperature  is  low,  the  rate  of  electrochemi¬ 
cal  reactions  in  the  reaction  layer  becomes  small,  and  the  current 
density  decreases  at  the  same  time.  For  the  positions  between 


x  =  0.01 5  m  and  0.02  m,  the  local  current  density  decreases  with  the 
increase  of  temperature.  This  is  mainly  due  to  the  lower  mole  frac¬ 
tion  of  hydrogen  and  oxygen  due  to  the  gradual  consumption  of 
hydrogen  and  oxygen  by  the  electrochemical  reactions. 

As  is  known  to  all,  the  average  current  density  directly  charac¬ 
terizes  the  overall  cell  performance.  When  the  operating  voltage  is 
kept  constant,  the  average  power  density  can  be  acquired  by  the 
production  of  average  current  density  and  operating  voltage.  Fig.  8 
shows  the  effects  of  emissivity  on  the  average  current  density  for 
a  co-flow  channel  arrangement  at  Vcen  =  0.75  V.  When  the  emissiv¬ 
ity  increases  from  0  to  0.2,  the  average  current  density  decreases 
rapidly,  and  with  the  further  increase  in  emissivity,  the  average 
current  density  reduces  gradually.  For  SOFCs,  the  surface  emissiv¬ 
ity  between  0.3  and  1  is  generally  adopted  by  various  researchers 
[28-31  ].  However,  from  Fig.  8  it  is  found  that  the  emissivity  varia¬ 
tion  between  0.3  and  1  has  a  small  influence  on  the  current  density. 
It  is  clear  that  the  numerical  results  rely  mainly  on  how  big  the 
emissivity  is. 

To  account  for  the  effects  of  emissivity  on  the  cell  electrochem¬ 
ical  performance,  a  relative  current  density  at  the  same  operating 
voltage  is  defined  as 

.,  Current  density  with  non-zero  emissivity 
—  Current  density  with  zero  emissivity 

with  which,  one  can  compare  the  magnitude  of  the  impact  of  the 
emissivity  on  the  average  current  density  at  different  operating 
voltage.  Fig.  9  gives  the  variation  of  the  relative  current  density 
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Fig.  5.  Local  temperature  distribution  in  the  plane  of  x  =  0.01 8  m  (close  to  the  outlet,  Co-flow,  unit:  K). 


with  emissivity  for  the  voltages  0.55  V,  0.75  V  and  0.85  V.  The 
trend  for  the  three  cases  is  the  same,  i.e.,  the  relative  current 
density  decreases  quickly  when  the  emissivity  is  between  0  and 
0.2,  and  after  that  it  varies  smoothly.  When  the  operating  voltage 
decreases,  the  average  current  density  and  the  heat  generation  will 
increase.  The  lower  the  operating  voltage,  the  greater  is  the  impact 


Fig.  6.  Effect  of  the  emissivity  on  the  axial  temperature  distribution  in  electrolyte 
(Vceii  =  0.75  V,  Co-flow). 


of  thermal  radiation.  However,  variation  of  the  current  density 
leads  to  a  variation  of  ohmic  heat  generation,  which  will  affect 
the  average  temperature  of  the  cell.  This  once  again  illustrates 
that  the  lower  the  operating  voltage,  the  greater  the  impact  of  the 
emissivity  on  the  average  temperature  of  cell.  This  is  confirmed 
by  Fig.  10.  In  addition,  the  net  radiative  heat  flux  is  affected  by 


Fig.  7.  Effect  of  the  emissivity  on  the  local  current  density  distribution  ( Vcen  =  0.75  V, 
Co-flow). 
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Fig.  8.  Effect  of  the  emissivity  on  the  average  current  density  (Vceii  =  0.75  V,  Co-flow). 


Fig.  9.  Effect  of  the  emissivity  on  the  relative  current  density  under  different  oper¬ 
ation  voltage  (Co-flow). 

the  cell  temperature.  Fig.  1 1  shows  the  effects  of  emissivity  on  the 
net  radiative  heat  flux.  It  is  found  that  the  net  radiative  heat  flux 
increases  rapidly  with  the  increase  of  emissivity  at  first.  However, 
the  increase  of  net  radiative  heat  flux  will  lead  to  a  decrease  of 


Fig.  10.  Effect  of  the  emissivity  on  the  cell  average  temperature  under  different 
operation  voltage  (Co-flow). 


Fig.  11.  Effect  of  the  emissivity  on  the  net  radiation  heat  flux  under  different  oper¬ 
ation  voltage  (Co-flow). 

the  cell  temperature,  and  in  turn  it  reduces  the  net  radiative  heat 
flux.  These  interactions  move  back  and  forth,  and  ultimately  a 
balance  between  temperature  and  net  radiation  heat  flux  might 
be  achieved,  which  leads  to  a  slight  variation  of  the  net  radiation 
heat  flux  and  cell  temperature  for  a  high  emissivity. 

3.2.  Influence  of  ambient  temperature  on  SOFC  performance 

Fig.  12  shows  an  interconnect  unit  adopted  in  the  experiment 
work  in  [32].  The  highlighted  region  is  the  single-cell  channel,  2 
and  3  are  the  inlet  and  outlet,  respectively.  The  complete  radia¬ 
tion  model  should  consider  the  radiative  heat  transfer  effects  of 
the  components  1  and  4.  However,  the  temperature  of  1  and  4  is 
difficult  to  obtain  from  the  single-cell  model,  and  must  be  solved 
through  the  cell-stack  model.  The  temperatures  of  1  and  4  might 
be  maintained  at  a  certain  value  if  a  suitable  thermal  management 
is  adopted  for  instance,  they  can  be  treated  as  the  environmental 
temperature.  In  Section  3.1 ,  the  ambient  temperature  is  set  equal  to 
inlet  temperature  of  the  fuel,  and  effects  of  different  ambient  tem¬ 
perature  on  cell  performance  will  be  discussed  in  this  Section.  Five 
cases  with  the  ambient  temperature  973  K,  1023  K,  1073  K,  1 123  K 
and  1173  K  are  investigated  separately. 


Fig.  12.  Interconnet  unit  in  Ref.  [32]. 
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Fig.  13.  Effect  of  the  ambient  temperature  on  the  average  net  radiation  heat  flux 
(e  =  0.5). 


The  effects  of  the  ambient  temperature  on  the  average  net  radi¬ 
ation  heat  flux  are  shown  in  Fig.  1 3.  When  the  ambient  temperature 
is  low,  the  average  net  radiation  heat  flux  is  positive,  which  means 
that  the  radiation  heat  transfer  rate  from  the  stack  to  the  ambient 
is  greater  than  that  from  ambient  to  the  stack.  With  the  increase 
of  ambient  temperature,  the  radiative  heat  transfer  flux  from  the 
ambient  to  the  cell  increases,  and  finally  it  exceeds  that  from  the 
stack  to  the  ambient,  which  leads  to  a  negative  net  radiation  heat 
flux.  In  addition,  at  a  constant  ambient  temperature,  a  lower  oper¬ 
ating  voltage  results  in  a  greater  net  radiation  heat  flux.  Fig.  14 
reveals  the  effects  of  ambient  temperature  on  the  average  cell  tem¬ 
perature.  It  is  found  that  the  average  cell  temperature  increases 
with  increasing  ambient  temperature.  When  the  operating  volt¬ 
age  is  kept  constant  (0.55  V,  0.75  V  and  0.85  V),  the  relative  average 
cell  temperature  for  Tamb  =  1 1 73  K  increases  about  2.83%,  2.80%  and 
2.78%,  or  32.1  K,  30.5  K  and  29.8  K,  respectively,  compared  to  those 
predicted  for  Tamb  =  973  K.  This  finding  reveals  that,  with  increasing 
ambient  temperature,  the  operating  voltage  has  a  small  influence 
on  the  average  cell  temperature.  The  effect  of  ambient  temperature 
on  the  average  current  density  is  illustrated  in  Fig.  15.  It  is  found 
that  increasing  the  ambient  temperature  has  a  small  effect  on  the 
average  current  density. 


Fig.  15.  Effect  of  ambient  temperature  on  the  average  current  density  ( e  =  0.5). 

3.3.  Influence  of  flow  pattern  on  SOFC  performance 

In  the  above  discussed  sections,  the  focus  is  on  the  performance 
of  SOFC  under  the  co-flow  condition.  Planar  stacks  can  be  charac¬ 
terized  according  to  arrangement  of  the  gas  flow:  Co-flow  (air  and 
fuel  flow  in  parallel  and  in  the  same  direction),  counter-flow  (in 
parallel  but  in  opposite  direction)  and  cross-flow  (perpendicular 
to  each  other).  There  is  a  great  difference  among  the  three  flow 
arrangements.  Only  the  counter-flow  and  the  co-flow  configura¬ 
tions  are  considered  in  this  study  because  numerous  meshes  are 
necessary  to  model  the  cross-flow  pattern  case  and  the  calculation 
time  becomes  very  long. 

At  the  same  operating  voltage,  the  effects  of  emissivity  on  the 
average  current  density  for  different  flow  arrangement  are  shown 
in  Fig.  16.  The  variation  of  the  current  density  for  the  counter-flow 
shows  the  same  trend  as  that  of  co-flow,  i.e.,  there  is  a  rapid  change 
in  the  current  density  when  the  emissivity  is  small,  and  the  varia¬ 
tion  tends  to  be  flat  when  the  emissivity  is  big.  Overall,  the  change 
in  counter-flow  is  smaller  than  that  in  co-flow  because  the  two 
flow  arrangements  have  different  local  temperature  distribution. 
Fig.  17  shows  the  local  temperature  distribution  in  the  air  channel. 
For  the  co-flow,  the  maximum  temperature  occurs  in  the  outlet 
region.  For  the  counter-flow,  a  constant  temperature  is  set  at  the 


Fig.  14.  Effect  of  the  ambient  temperature  on  the  average  temperature  of  cell 
(£  =  0.5). 


Fig.  16.  Effect  of  the  emissivity  on  the  average  current  density  under  different  flow 
arrangement  (Vcen  =  0.75  V). 
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inlets  of  air  and  fuel,  which  leads  to  a  maximum  temperature  in 
the  middle  part  of  the  cell.  Because  the  view  factors  between  the 
ambient  and  surfaces  near  the  inlet  and  the  outlet  are  much  bigger, 
the  corresponding  radiative  heat  transfer  is  also  greater.  On  the 
other  hand,  the  radiative  heat  transfer  between  the  ambient  and 
the  surface  near  the  middle  part  of  the  cell  is  lower.  Considering  the 
maximum  temperatures  for  the  co-flow  and  the  counter-flow,  it  is 
clearly  found  that  the  emissivity  has  less  influence  on  the  current 
density  for  counter-flow  arrangement.  The  effects  of  emissivity  on 
the  net  radiation  heat  flux  for  different  flow  arrangement  are  shown 
in  Fig.  18.  The  emissivity  also  has  less  influence  on  the  average  net 
radiation  heat  flux  for  counter-flow  arrangement. 

Fig.  19  shows  the  effect  of  the  ambient  temperature  on  the 
current  density.  At  the  same  operating  voltage,  the  ambient  tem¬ 
perature  has  less  influence  on  the  current  density  for  the  counter¬ 
flow,  compared  to  the  co-flow  arrangement. 


Fig.  17.  Local  temperature  distribution  in  air  channel  under  different  flow  arrange¬ 
ment  (Vceii  =  0.75  V,  e  =  0.5). 


Emissivity  e 


4.  Conclusions 

To  investigate  radiative  heat  transfer,  a  3D  comprehensive 
model  has  been  developed  for  calculation  of  the  performance  char¬ 
acteristics  of  planar  solid  oxide  fuel  cells.  The  effects  of  emissivity, 
ambient  temperature  and  flow  arrangement  are  investigated  in 
details.  The  major  findings  are  as  follows. 

(1)  With  the  increase  of  emissivity,  the  current  density  decreases 
quickly  at  small  s  but  only  slightly  at  big  s ,  particularly  at  low 
operating  voltage. 

(2)  For  the  typical  materials  for  SOFC  the  emissivity  is  between  0.3 
and  1,  and  it  is  found  that  the  emissivity  has  a  small  influence 
on  the  simulation  results  when  thermal  radiation  is  considered. 

(3)  The  effects  of  ambient  temperature  on  the  performance  of 
SOFC  are  relatively  small.  With  increasing  ambient  tempera¬ 
ture,  the  average  temperature  and  current  density  of  the  cell 
will  increase,  but  the  variation  is  small.  In  addition,  the  relative 
current  density  increases  faster  than  the  average  temperature. 

(4)  For  different  flow  arrangements,  the  radiation  for  co-flow  has  a 
bigger  influence  on  the  performance  of  SOFC  than  for  counter¬ 
flow. 


Fig.  18.  Effect  of  the  emissivity  on  the  net  radiation  heat  flux  under  different  flow 
arrangement  ( Vcen  =  0.75  V). 


Fig.  19.  Effect  of  the  ambient  temperature  on  the  average  current  density  under 
different  flow  arrangement  (Vcen  =  0.75  V,  e  =  0.5). 
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